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ABSTRACT 



We study semianalytically the gravitational effects on neutrino pair 
annihilation near the neutrinosphere and around the thin accretion disk. For 
the disk case, we assume that the accretion disk is isothermal and that the 
gravitational field is dominated by the Schwarzschild black hole. General 
relativistic effects are studied only near the rotation axis. The energy deposition 
rate is enhanced by the effect of orbital bending toward the center. However, the 
effects of the redshift and gravitational trapping of the deposited energy reduce 
the effective energy of the gamma ray bursts' source. Although each effect is 
substantial, the effects partly cancel one another. As a result, the gravitational 
effects do not substantially change the energy deposition rate for either the 
spherical symmetric case or the disk case. 

Subject headings: accretion, accretion disks — black hole physics — gamma rays: 
bursts 
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1. INTRODUCTION 

The relativistic fireball (Shemi k Piran 1990; Rees k Meszaros 1992; 
Meszaros k Rees 1993; Sari k Piran 1995; Sari, Narayan k Piran 1996) is one 
of the most promising models of gamma ray bursts (GRBs). However, even if 
the fraction of the baryon rest energy is only 10 -3 in the fireball, the relativistic 
bulk flow, which is indispensable to GRBs, cannot be realized. Notwithstanding 
the very high energy phenomenon (10 52 ergs), the baryon density in the fireball 
must be extremely small. This is the famous baryon contamination problem 
and still remains unsolved. Thus the central engine of GRBs is still beyond 
deep mist. The source of GRBs may be one super massive (failed) supernovae 
(Woosley 1993; Paczyhski 1998) or may be a merger of two neutron stars or of a 
neutron star and a black hole (e.g. Eichler et al. 1989; Narayan, Paczyhski k 
Piran 1992; Meszaros k Rees 1992a; Katz 1997; Ruffert k Janka 1998, 1999). 

In these compact high energy objects, the neutrino-antineutrino annihilation 
into electrons and positrons (hereafter neutrino pair annihilation) is a possible 
and important candidate to explain the energy source of GRBs (Paczyhski 
1990; Meszaros k Rees 1992b; Janka k Ruffert 1996; Ruffert et al. 1997; 
Ruffert k Janka 1998, 1999). Motivated by the delayed explosion of Type II 
supernovae, the energy deposition rate due to the neutrino pair annihilation 
above the neutrinosphere has been calculated (Goodman, Dar k Nussinov 1987; 
Cooperstein, Van Den Horn k Baron 1987; Berezinsky k Prilutsky 1987). The 
energy deposition rate is proportional to r~ 8 (r is the distance from the center 
of the neutrinosphere) for a large r, and almost all deposition occurs near the 
neutrinosphere. As they themselves noted in their paper, Goodman et al. (1987) 
neglected the general relativistic effects on the energy deposition rate, which 
may change their numerical value seriously. In simulations of the neutrino pair 
annihilation rate, it is very important to confirm whether or not the energy 
deposition rate is altered or not by the gravitational effects. In the recent study, 
Salmonson k Wilson (1999) concluded that the energy deposition rate in Type 
II supernovae is enhanced about 4 times as a result of the gravitational effects. 
We must check whether or not their results can be applied to the central engine 
of GRBs. 

One of the most probable candidates for the central engine of GRBs is the 
accretion disk around a black hole (Woosley 1993; Popham, Woosley k Fryer 
1999; MacFadyen k Woosley 1999; Ruffert k Janka 1999). The system of an 
accretion disk and a black hole may be formed by the merging of two neutron 
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stars, the merging of a black hole and a neutron star, or the failed supernovae. 
In general, the baryon density has the lowest value along the rotation axis just 
above the black hole (e.g. see Ruffert & Janka 1999). This region might be a key 
to resolving the baryon contamination problem. The hot accretion disk emits 
neutrinos and antineutrinos. The energy deposited in the lowest density region 
is a candidate for the central engine of GRBs. Using hydrodynamic simulations, 
Ruffert & Janka (1999) showed that the neutrino pair annihilation deposits 
energy in the vicinity of the torus at a rate of (3 — 5) x 10 50 ergs s^ 1 . They 
concluded that the gravitational effect on the energy deposition rate around the 
accretion disk is small. We must supplement their results from the analytical 
side. 

In various arguments on the energy deposition in the central engine of GRBs, 
the order estimation of the deposited energy is sufficient, at least at present. In 
this article, based on simple models, we study semianalytically the gravitational 
effects on the energy deposition rate for two cases. In one case neutrinos are 
emitted spherically symmetrically. In the other case the hot accretion disk 
emits neutrinos. We have derived the gravitational effects on the former case 
independently of Salmonson & Wilson (1999). Some differences of our work 
from Salmonson & Wilson in the formulation, the interpretation of the energy 
deposition, and the additional factor are mentioned. As for the disk case, we 
assume that the accretion disk is isothermal and that the gravitational field is 
dominated solely by the central Schwarzschild black hole. These assumptions 
enable us to treat the energy deposition around the disk semianalytically. Thus, 
in both two cases gravitation is described by the Schwarzschild metric, and the 
essential differences between the two cases come from the shape of the neutrino 
emitters. 

The gravitational effects consist of three factors: they are the bending of 
neutrino trajectories, the gravitational redshift, and the trapping of deposited 
energy into the central gravitational source. We show that the energy deposition 
rate is indeed enhanced rather crucially by the effect of neutrino bending. 
However, it is also shown that the gravitational redshift and the trapping of the 
deposited energy reduce this enhancement. As a result, the gravitational effects 
do not substantially change the energy deposition rate for either the spherical 
symmetric case or the disk case. 

This paper is organized as follows. In section two we investigate neutrino pair 
annihilation near the neutrinosphere. The same process around the accretion 
disk is discussed in section three. The last section is devoted to conclusions. 
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2. NEUTRINO PAIR ANNIHILATION NEAR THE 
NEUTRINOSPHERE 

In this section we study the general relativistic effects on neutrino pair 
annihilation near the neutrinosphere. This study has been already done by 
Salmonson & Wilson (1999). Using another method, we formulate the same 
problem independently of the work of Salmonson & Wilson. Some alterations 
in the interpretation of the energy deposition in Salmonson & Wilson are 
mentioned. 

The number of reaction, v + v — > e + + e~, per unit volume per unit time 
(Goodman, Dar & Nussinov 1987) is written as 

^tdV = I J f»( p ^ r W*(P^ r "> a \ v »~ v »\ d3 P» d3 P»- W 

Here /„ is the number density of neutrinos (antineutrinos) in phase space, 
v u (vp) is the velocity of neutrinos (antineutrinos), and a is the rest-frame 
cross section. The left handside of equation (P is Lorentz invariant, since both 
the numerator, dN, and denominator, dtdV = y/—gd 4 x, are Lorentz invariant. 
Since f v and d 3 p u /e u (where e u is the proper energy of neutrinos) of the right 
handside are also Lorentz invariant, e u e 9 \v v — v p \a should be Lorentz invariant. 
The latter is written in a manifest Lorentz- invariant form as ac 3 (p u ■ p„), where 
(Pv ' Pv) is the inner product of the 4-momenta. The standard model predicts 
that the cross section is expressed as 

a = 2c 2 KG 2 F (pvPp), (2) 
where the dimensionless parameter K is written as 



K{y^y^) = K{v T v T ) 



1 - 4 sin 2 #w + 8 sin 4 ^ 



6tt 



„, l+4sin 2 #w + 8sin 4 # w 

K(u e v e ) = . (3) 

Here the Fermi constant Gp = 5.29 x 10~ 44 cm 2 MeV -2 and the Weinberg angle 
sin 2 fl w = 0.23. 

Let us incorporate the effects of gravitational force due to the neutron 
star or black hole on the neutrino pair annihilation rate. We assume that the 
gravitational field is described by the Schwarzschild metric: 

ds 2 = g ij dx i dx j = (l-^ c 2 dt 2 - ^-^dr 2 - r 2 (d6 2 - sin 2 6d^ 2 ) , (4) 
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where r g = 2GM/c 2 is the Schwarzschild radius. In this field the eikonal for a 
massless particle (Landau & Lifshitz 1979) is written as 

■ip = —uj t + Lip + ip r (r), (5) 

where u> and L are constants. Vv( r ) satisfies the equation 

dVv(r) ujI f r g y 2 L 2 1 
<9r \ c 2 \ r ) r 2 l- r f K ' 

From equation (|5|), we can obtain the momentum of a neutrino by p, = ftj^- 

Let us consider a neutrino and an antineutrino moving on the same surface, 
9 = 7r/2. In this case, the inner product of the two particles is written by 

ipvPs) = 9 l3 PviPvj (7) 

= £ -^f(l y/l fMVl 

P^Pp ^ _ rp 



where 



The proper energy of the neutrino has been written as 



1- — — , 1-*' 



(10) 



where Eq u is the energy observed at infinity. Thus the proper energy is redshifted, 
as is well known. If we define an angle 9 V as 



sin0„ = ^Jl-^, (11) 
f y f 

equation (§) becomes a simple and natural form, 

(pvPv) = £ -^f(l-cos(9v-9u)). (12) 

The angle 9 V (Op) represents the angle between p v (p p ) and the position 
vector r (see Figure 1). We assume that the neutrinosphere emits neutrinos 
and antineutrinos isotropically. Then we can write the number densities as 
fv(p u ,f)d 3 p u = n{e u )e 2 de u dVL. Because p u is constant along a neutrino ray, the 
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maximum angle, 8m, is obtained by substituting tt/2 for Q v at the radius of the 
neutrinosphere, R u , in equation (|TT|). Thus we obtain 



sin 6^m = — 

r 



\ 



1 - ^ 

Y^±- (13) 



The effect of the orbital bending is apparent in this equation. Until now we 
have discussed the maximum angle on the surface of 9 = it/2. In general cases, 
the angles between p u and r or the inner product, {p u ■ p p ), are expressed by the 
two angles, 6 V and <p u . From the symmetry, the behaviour of 9m is obviously the 



same as described in equation (|i~3l) , and <p u varies from to 2ir. 

Using the effective temperature of the neutrinosphere, T eff = To/y^oo with a 
constant T , we can write the density 

n{e v ) = j^-r j^-^ , (14) 

where g u is a statistical factor (g v = 1 for a neutrino). e u /(kT e ff) is constant 
along a neutrino ray, since the redshift is cancelled out. Thus n(e u ) is conserved 
along a neutrino ray in accordance with Liouville's theorem in curved spacetime 
(Misner, Thorne & Wheeler 1975). From the above formulation, one can find 



dN(r) 



2cKGlF(r) J J de v de D n{e u )n{e D )elel, (15) 



dtdV 

where the dimensionless factor F(r) is written by 

r&M pOm r 27r r 27r 

F(r) = / d6 u sin Q v I d9 D sin 9 P / dip v I d(p t 
Jo Jo Jo Jo 



x (1 — sin d v sin Q v cos (ip u — (p B ) — cos 6 V cos 6 P ) 2 (16) 
■(X 2 + 4X + 5), (17) 



27^(1 -X)* v , 



where 



X 



\ 



R \ 2 1 - ^ 
1"(-)t-^- (18) 



r 



' 9 



In our assumption, the energy deposited by the neutrino pair annihilation is 
propagated outward as a fireball or a shock wave, and observed as a GRB by a 
distant observer. Thus the energy we need to calculate is e 0l/ , not the proper 
energy e v . In this case the energy deposition rate is obtained by putting a factor 
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(e 0u + e p) i n the integrand in equation 

dE (r) 2cKG\ r<*> ; x 



dtdV (i_^yJo Jo 

xn(e v)n(e 0D )el v el P (e 0v + e 0P )F(r) (19) 



2 ^^) K wf> (2°) 



(1 - ?) 



The integrals for Eq u arid £op should be defined in the range in which the 
total energy produced by the pair annihilation is larger than the mass of 
created electrons, and smaller than the masses of weak bosons. Here we have 
approximated the integrals as expressed in equation (|T9|) in the same manner 
as Salmonson & Wilson (1999) did. This is because the cross section decreases 
with the energy of neutrinos, and the number of neutrinos whose energy is larger 
than the masses of weak bosons is also very small in our assumption (kT e ff is of 
the order of several MeV). The factor (1 — r g j R„) 9 ^ 2 /(l — r 3 /r) 4 represents the 
effect of the gravitational redshift, and F(r) includes the effect of the orbital 
bending. 

As is understood from the Lorentz invariant, dtdV = ^/—gd 4 x, if we integrate 
equation (p0|) over proper volume, dV = ^ —g Tr g e Qg w drdQdtp, we can obtain 
the total energy deposition per unit proper time, dr = i/goodf. It is natural 
to evaluate the energy deposition rate by the world time dt for a distant 
observer. We integrate over the volume, dV = ^jgeeg^drdOdtp. Thus the energy 
deposition per unit world time is expressed as 



dE Zlit* KGjgj 9 f r g \2 r<*> 2 F(r) 

-IT = — a5) n^r {kTcs) V - r-J L dr47Tr (TT^7 W ' { ) 

where we have put a factor, 

C(r) -Ui (I -:h 122) 




in the integrand. This is the escape probability of the deposited energy at r 
from the gravitational attraction (Chandrasekhar 1983; Shapiro & Teukolsky 
1983; Ruffert & Janka 1999). The electrons, positrons and photons which are 
captured by the gravitational attraction cannot contribute to the energy source 
of GRBs. Apart from C(r), the radial profile in the integrand of equation (pip is 
different from those of Salmonson & Wilson (1999), since Salmonson & Wilson 
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calculated the proper energy deposition per unit proper time. Of course, the 
results of Salmonson & Wilson are not mistakes for the estimate of the energy 
deposition rate in supernovae. For the source of GRBs, however, equation ( pTj ) 
is adequate. 

Let us investigate the effects of the redshift, orbital bending and gravitational 
capture. We integrate equation fl2~If ) and obtain the energy deposition rate for 

v e as 

^ = 1.27xl0 42 f^Vf4^V/ ergss-\ (23) 
dt \lMeVJ viow J B K J 

where the dimensionless factor / expresses the effects of the general relativity 
(/ = 1 when we neglect the gravitation). The energy deposition rates for 
and u T are 0.64 times equation ([23|). We numerically estimate / including the 
effects of the redshift only, the orbital bending only, or both redshift and orbital 
bending. Last, the total effects of the redshift, orbital bending, and gravitational 
capture are calculated. The results are listed in Table 1. As Table 1 or equation 
( pID indicates, the effect of the redshift reduces the energy deposition rate, and 
the effect of the orbital bending increases it. Although each effect, that of the 
redshift and that of orbital bending, is substantial, the effects partly cancel each 
other. As a result, the order of the energy deposition rate for the most probable 
case, R u /r g = 2.5, is not altered. When we neglect the general relativistic 
effects, the energy deposition rate increases by 1.3 times as R v becomes 10% 
larger, and also increases by 2.4 times as the temperature becomes 10% higher. 
Therefore, the gravitational effects are not so large in comparison with the 
errors due to the uncertainties of R v or T e g, and are overwhelmed by them. The 
effect of the gravitational capture becomes important as R u /r g decreases. As is 
plotted in Figure 2, in the cases of both the presence and absence of gravitation, 
the energy deposition mainly occurs near the neutrinosphere. 

Salmonson Sz Wilson (1999) concluded that the effects of gravity enhance the 
energy deposition rate up to a factor of more than 4 for R u < 2.5r g . However, 
our results show that the gravitational effects reduce the energy deposition 
rate. This discrepancy survives even if we omit the escape factor C(r). The 
proper energy deposition per unit proper time is enhanced by both the effects 
of the redshift and that of orbital bending. Additionally, Salmonson & Wilson 
expressed the general relativistic effects with the fixed neutrino luminosity at 
infinity Loo, whereas we have done so with the local physical quantity T e g. 
Therefore, an additional factor coming from the redshift of the luminosity 
(L(R V ) oc T e 4 ff , Lqo = (1 — r g /Ru)L(R u )) enhances the energy deposition rate in 
the work of Salmonson & Wilson. However, the quantity L^o of GRBs is not 



-10- 



directly observable at present. It is more natural to study the effects for the 
given local parameters, T e g or L(R U ), which is restricted or provided by models 
of the central engine. 



3. NEUTRINO PAIR ANNIHILATION AROUND THE 

ACCRETION DISK 

In this section we investigate the energy deposition rate around the accretion 
disk. In order to simplify our formulation, we assume that the accretion disk 
is isothermal and that the gravitational field is dominated by the central 
Schwarzschild black hole. We neglect the rotation of the black hole. The 
accretion disk is assumed to be thin, and its self-gravitational effects are 
neglected. Of course, these idealizations may be far from the case of the realistic 
accretion disk. However, we consider that this simple method is sufficient for 
qualitatively studying the gravitational effects on the energy deposition rate. In 
this case the equation of the energy deposition rate is the same as equation (f2~0|) 
provided that F(r) is replaced by F(r,0) (it will be given below). The effect of 
the gravitational redshift can be easily incorporated, whereas the formulation of 
the neutrino bending is difficult to do because the accretion disk emits neutrinos 
anisotropically. 

First, we calculate the dimensionless factor F(r, 9) without the effect of 
gravity. The accretion disk is placed on the equatorial plane, 9 = n/2. The 
black hole is at the origin, and we consider a point P = (r, 9, 0) where pair 
annihilations occur (see Figure 3). A neutrino is emitted from an arbitrary 
point on the disk S = (R, 7r/2, (p), where R is limited in the range from R[ n to 
-R out . The neutrino emitted from 5" travels straight and arrives at the point P. 
Let us denote the angle components of the vector joining S and P by (9 U , <p v ). 
They are given by 

r cos 9 , . 

V ' + R — 2rR sm 9 cos <p 

—R sin if 

smip u = (25) 
J r 2 sin 2 9 + R? — 2rR sin 9 cos cp 

Thus 9 U and <p v are functions of R and (p for fixed r and 9. The Jacobian 
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J = d{6 u1 ip u )/d(R,ip) is 

rR cos 9 , , 

J = . (26) 

'r 2 sin 9 + R 2 — 2rR sin 9 cos ip(r 2 + R 2 — 2rR sin 9 cos <£>) 
Consequently, we obtain F(r,9) as 

rRout rRout r^ r^n 

F(r, 9) = dR dR' / d<p / ci^'JJ' 

JR in JR in JO JO 

x sin 9 V sin P (1 — sin # v sin 6* p cos — (p p ) — cos 6^ cos 9 P ) 2 . (27) 



In equation ([270 we adopt i? in = 3r g , the innermost stable orbit, and 
-Rout = 10r 9 as Woosely (1993) assumed. F(r, 9) derived from the numerical 
integral of equation ( p7|) is plotted in Figure 4(a) and (b). As is shown in these 
figures, the energy deposition rate is maximized in the vicinity of the accretion 
disk, where F(r, 9) ~ 30 — 33. The simulation of a neutron star merger by 
Ruffert & Janka (1999) showed that the Paczyhski-Wiita potential (Paczyhski & 
Wiita 1980), which mimics the effects of the general relativity, gives a relatively 
more transparent disk for neutrinos than that given by the Newtonian potential. 
The profile of the energy deposition rate in the Paczyhski-Wiita potential is 
similar to our analytical one depicted in Figure 4(a), which shows that the 
rate takes its maximum value on the surface of the disk. On the other hand, 
the simulated deposition rate in the Newtonian potential is maximized near 
the rotation axis. Let us calculate the energy deposition rate near the rotation 
axis, where is the lowest baryon density region. Thus we calculate in the region 
9 < 7r/4 and obtain 



where the dimensionless quantity G shows the relative contributions from 
various regions in the absence of gravitation. G is normalized to unity when 
we integrate over the volume for 9 < 7r/4 and r = 2r g — 10r g . The values of G 
in the other regions are summarized in Table 2, from which we can obtain the 
energy deposition rate in the respective region. We neglect the energy deposited 
inside r = 2r g , since the baryon density in this region is very high and the 
energy contribution is small for the small volume and deposition rate. In the 
case of the spherical emitter in section 2, the deposition rate decreases as r~ 8 . 
On the other hand, around the accretion disk, as is seen from Table 2, there 
remains a marginal deposition rate even at regions relatively distant from the 
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center. Of course the deposition rate per unit volume at distant positions is 
small. However, large volume results in a non-negligible contribution at the 
regions distant from the center. 

Untill now we have neglected the gravitational effects. It is easy to 
incorporate the effects of the redshift and trapping by the central gravitational 
source in the preceding arguments of this section. However, the bending effect is 
difficult to treat, unlike the case of the neutrinosphere, since the accretion disk 
emits neutrinos anisotropically. Thus we are forced to make some approximation. 
As is shown in Figure 4(b), the 9 dependence of F(r,9) is weak for small 9. 
We may set F(r,8) ~ F(r, 0) for 9 < 7r/4. In the absence of gravitation if 
we adopt this approximation in the region 9 < rr/4 and r = 2r g — 10r g , we 
obtain G = 0.81. The exact value of G is unity, and this approximation is not 
necessarily satisfactory. However, this approximation may be sufficient for the 
order estimate of the gravitational effects. 

We can obtain F(r, 0) including the effect of orbital bending with comparative 
ease, since the geometry of this case maintain the symmetry. A neutrino is 
emitted from the disk at R and 9 = n/2, and it arrives at a point at r and 9 = 0. 
The nearest distance, ro, from the origin to the orbit of neutrinos (Landau & 
Lifshitz 1979) is numerically obtained from 

r dr' . . 

*/ 2 = L 1 o (29) 



c 
r 



Here, in the case in which a neutrino passes through r until it arrives at a point 
at 9 = 0, the integration for r' is performed from ro to R and r. When the 
distance from the origin to the neutrino varies monotonically, the integration is 
performed from the smaller to the larger of r and R. We can get 9 V at 9 = 
numerically from r and the following equation; 



sine',, = — 



r 



\ 



ro 



(30) 



The constant, ro, or 9 U , is a function of r and R in this case. As is easily 
understood, a neutrino coming from R[ n forms 9 m , the minimum value of 9 V , 
at 9 = and that from R out forms the maximum value of 9, 9m- Integrating 
equation fllED from 9 m to 6*m, we obtain F(r, 0) involving gravitational effects. 
In Figure 5 we plot F(r, 0) for both the case when the bending is taken into 
consideration and the case when it is not. In comparison with the spherical case 
in Figure 2, the deposited energy at distant regions in the disk case is marginally 
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substantial. In the presence of bending, the peak of the energy deposition rate 
is shifted to a little bit larger r and the value of the rate at the peak is about 
twice as larger as values obtained in the absence of bending. 

Although the 9 dependence of the gravitational effects may not necessarily 
be small, unlike F(r, 9) in the absence of gravitation, we assume it is small 
here. Using F(r, 0) with the bending effect, we calculate the energy deposition 
rate in the range 9 < n/A and r = 2r g — I0r g . Table 3 lists the values of the 
factor / that shows the gravitational effects. The ^-dependence is neglected 
in our calculation except for the case involving the redshift only. Thus / is 
normalized to unity when F(r, 9) (in the case involving the effect of the redshift 
only) or F(r, 0) (in the other cases) is integrated over r and 9 in the absence of 
gravitation. As is easily seen, the gravitational effects cancel one another out. 
This is analogous to the neutrino sphere case in the previous section. This result 
strongly supports that of Ruffert & Janka (1999). They treated the system of 
the accretion torus and a black hole unlike our system of the disk and a black 
hole. Using an approximation similar to ours, they analytically calculated the 
energy deposition rate due to the neutrino pair annihilation. Their result is that 
the gravitational effects reduce the deposition rate by a factor of 10 — 30%. It 
agrees well with our result. 

In order to circumvent the baryon contamination problem, the energy 
fraction of baryonic matter in the fireball must be less than about 10~ 5 (Shemi 
& Piran 1990). If we adopt the duration time of the neutrino radiation to 
be tdur = 0.1s and T c g = lOMeV, the highest mean mass densities p inside 
9 = — 7r/3 to resolve the above problem are 10 6 g/cm 3 for r = 2r g — 5r g , 
10 5 g/cm 3 for r = 5r g — 10r g and 10 4 g/cm 3 for r = 10r 9 — 20r g . Since some 
fraction of energy really escapes from the considered regions during the finite 
duration time, the above restrictions may become more stringent. 

4. CONCLUSIONS 

In this article we have investigated semianalytically the neutrino pair 
annihilation near the neutrinosphere and around the thin accretion disk 
assuming that the gravitational sources in both cases are described by the 
Schwarzschild metric. The accretion disk has been assumed to be a blackbody 
and isothermal. These assumptions enable us to treat these two cases based in 
an almost unified fashion, which also clarifies the physical differences between 
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these two cases. We have studied the general relativistic effects only near the 
rotation axis, because that region is especially of interest to the source of GRBs 
and estimating the effect of orbital bending for large 6 is difficult. 

The general relativistic effects as a whole do not enhance the neutrino 
energy deposition rate in either case. The energy deposition rate is enhanced 
by the effect of orbital bending toward the center. However, the enhancement 
is cancelled out by the effects of the redshift and capture by the gravitational 
attraction. Consequently, numerical simulations of the neutrino energy 
deposition rate in various models can correctly estimate the order of the rate 
without considering the gravitational effects, since it is supposed that the 
thickness, shape, or temperature distribution of the disk or sphere does not 
greatly affect the gravitational effects themselves. Taking into account also the 
results of Ruffert & Janka (1999), the conclusions mentioned above are strongly 
suggested to be valid in the following geometrical forms of the neutrino source: 
sphere, thin disk and torus. We have also shown that the neutrinos emitted 
from the disk can deposit energy at more distant regions than the neutrinos 
emitted from the sphere. The importance in this article resides in the qualitative 
properties of the general relativistic effects. The quantitative calculations in 
this paper are not so important, and should be investigated on the basis of more 
sophisticated models and simulations. 
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Figure captions 

Fig. 1. — Maximum angle 9m at a distance r from the center formed by neutrinos 
emitted from the surface of the neutrinosphere. The dashed line shows the orbit 
forming the maximum angle in the absence of gravitation. 

Fig. 2. — The plot of F(r) against r for R v = 2.5r g when neutrinos are emitted 
isotropically. The solid (dashed) line shows F(r) for the case in the presence 
(absence) of gravitation. The origin is at r = R v . 

Fig. 3. — Geometry of neutrino's path. A neutrino is emitted from a point 
S = (R, 7r/2, <f) on the disk and arrives at P = (r, 9, 0). The angular components 
of the vector joining S and P are {Q v ,(p v ). R (radius on the disk) is limited to 
the range from R in to R out - 

Fig. 4. — Contour plot of F(r, 9) where neutrinos are emitted from the disk. Here 
gravitation has been neglected, (a) Plot in Cartesian coordinates (r sin 9, r cos 9). 
(b) Plot in polar coordinates, (r,9). The contours are plotted at unit intervals 
except for the line of F(r,9) = 0.1; 9 = and 9 = n/2 correspond to points on 
the rotation axis and on the disk, respectively. 

Fig. 5. — Diagram of F(r, 0) vs. r where neutrinos are emitted from the disk. 
The solid (dashed) line shows F(r, 0) for the case in which we consider (neglect) 
the effect of gravitational bending. 
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Redshift Only 


Bending Only 


Redshift and Bending 


Whole 
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0.32 


4.7 


0.97 


0.57 


2.5 


0.60 
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0.87 


0.73 
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0.81 
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0.93 
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Table 1. The dimensionless factor / represents the general relativistic effects 



(see eq. [23]) when neutrinos are emitted isotropically from the neutrinosphere; 
/ is normalized to unity in the absence of gravitation. The column headings 
"Redshift Only" and so on indicate the incorporated effects of gravitation; 
"Whole" over the last column means that we incorporate all gravitational 
effects, redshift, bending, and trapping. The energy deposition rate is enhanced 
by orbital bending and reduced by the redshift and trapping. 



9 2r g — 5r g 


5r g — 10r g 


10r 9 - 20r 9 


- tt/4 0.35 


0.65 


0.22 


7r/4-7r/3 0.32 


0.77 


0.17 



Table 2. The dimensionless factor G. It represents the fraction of the energy 
deposition rate for each region (see eq. [28]) when neutrinos are emitted from 



the disk. G is normalized to unity for the region surrounded by r = 2r g — 10r g 
and 9 = — it/ 4. Here we neglect the effects of gravitation. 



range / 

redshift only bending only redshift and bending whole 

2r g -hr g 0.66 2.6 1.6 1.4 

5r g -l0r g 0.31 2.5 0.78 0.75 

Table 3. The dimensionless factor / for 9 = — tt/4 (see equation ||28|| ). 
Neutrinos are emitted from the disk; / is normalized to unity when F(r, 9) (in 
the case involving the effect of the redshift only) or F(r, 0) (in the other cases) 
is integrated over r and 9 in the absence of gravitation. 
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